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Quantum Monte Carlo calculations of electromagnetic moments and transitions are reported for 
A < 9 nuclei. The realistic Argonne vis two-nucleon and Illinois-7 three-nucleon potentials are used 
to generate the nuclear wave functions. Contributions of two-body meson-exchange current (MEC) 
operators are included for magnetic moments and Ml transitions. The MEC operators have been 
derived in both a standard nuclear physics approach and a chiral effective field theory formulation 
with pions and nucleons including up to one-loop corrections. The two-body MEC contributions 
provide significant corrections and lead to very good agreement with experiment. Their effect is 
particularly pronounced in the A — 9, T = 3/2 systems, in which they provide up to ~ 20% (~ 40%) 
of the total predicted value for the 9 Li ( 9 C) magnetic moment. 
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I. INTRODUCTION 

Quantum Monte Carlo (QMC) calculations of elec- 
troweak transitions in A — 6,7 nuclei were reported 
in Ref. [l| and corrections for the magnetic moments 
(m.m.'s) and Ml transitions from two-body meson- 
exchange current (MEC) operators were given in Ref. . 
The QMC method is a two-step process, with an ini- 
tial variational Monte Carlo (VMC) calculation to find 
a good trial function, followed by a Green's function 
Monte Carlo (GFMC) calculation to refine the solution. 
When used with the Argonne vis two-nucleon 0] and 
Illinois-2 three-nucleon [4| potentials, the final GFMC 
results reproduce the ground- and excited-state energies 
for A < 10 nuclei very well [B-Q. 

In the present paper, we extend these calculations 
to A = 8,9 nuclei using the improved Illinois-7 three- 
nucleon potential @. The electromagnetic (EM) cur- 
rent operator includes, in addition to the standard one- 
body convection and spin-magnetization terms for in- 
dividual protons and neutrons, a two-body MEC com- 
ponent. The latter is constructed within two distinct 
frameworks, namely the same standard nuclear physics 
approach (SNPA) illustrated in Refs. and the chi- 

ral effective field theory (xEFT) formulation of Refs. fllT - 

We report energies, radii, magnetic and quadrupolc 
moments, and a number of Ml and E2 transitions. The 
MEC contributions can make significant corrections to 
the m.m.'s and Ml transitions, and we find general agree- 
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ment between the two formulations and with experiment. 
However the ;\EFT formulation provides better agree- 
ment for the calculated m.m.'s, for which both MEC 
models are tested. The Ml transitions are calculated 
only with the %EFT MEC operators, showing improved 
agreement with experiment in all cases. 

A brief review of the QMC calculational method is 
given in Sec. [TXJ The EM current operator is discussed 
in Sec. HID Results and conclusions are given in Sees. IIVI 
and ED 



II. QUANTUM MONTE CARLO METHOD 

We seek accurate solutions of the many-nuclcon 
Schrodinger equation 

H^>{r;T,T z ) = E*{r-T,T Z ) , (1) 

where ^(J 71 ;T,T Z ) is a nuclear wave function with spe- 
cific spin-parity J*', isospin T, and charge state T z . The 
Hamiltonian used here has the form 

ff = E^+X>«+ E v v>> > ( 2 ) 

i i<j i<j<k 

where Ki is the non-relativistic kinetic energy and 
and Vijk are respectively the Argonne vis (AV18) [|[ and 
Illinois-7 (IL7) @ potentials. 

The VMC trial function 4v( J v ; T, T z ) for a given nu- 
cleus is constructed from products of two- and three-body 
correlation operators acting on an antisymmetric single- 
particle state of the appropriate quantum numbers. The 
correlation operators are designed to reflect the influence 
of the interactions at short distances, while appropriate 
boundary conditions are imposed at long range Hid- 
The ( J ff ; T, T z ) has embedded variational parameters 
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that are adjusted to minimize the expectation value 



(3) 



which is evaluated by Metropolis Monte Carlo integra- 
tion . Here Eq is the exact lowest eigenvalue of H for 
the specified quantum numbers. A good variational trial 
function can be constructed with 
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where the S is a symmetrization operator. The Jas- 
trow wave function tyj is fully antisymmetric and has 
the (J^;T,T Z ) quantum numbers of the state of in- 
terest, while Uij and Uf^ 1 arc the two- and three- 
body correlation operators. Although we construct the 
\&y( J 71 '; T, T z ) to be an eigenstate of the isospin T, we 
allow isobaric analog states with different T z to have dif- 
ferent wave functions, reflecting primarily the difference 
in Coulomb contributions, but also additional charge- 
symmetry-brcaking parts of the AV18 interaction. 

The GFMC method 0, [H| improves on the VMC 
wave functions by acting on Sty with the operator 
exp [— (H — Eq) t\. In practice, a simplified version H' 
of the Hamiltonian H is used in the operator, which in- 
cludes the isoscalar part of the kinetic energy, a charge- 
independent eight-operator projection of AV18 called 
AV8', a strength-adjusted version of the three-nucleon 
potential IL7' (adjusted so that (H') ~ (H)), and 
an isoscalar Coulomb term that integrates to the total 
charge of the given nucleus. More detail can be found in 
Refs. [lEIli- 

The operator is applied in small slices of imaginary 
time r to produce a propagated wave function: 
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Obviously \P(r = 0) = v&y and *(t -> oo) = * . The 
algorithm for propagation produces samples of the wave 
function ^(r) but does not provide gradient information. 
Therefore, quantities of interest are evaluated in terms of 
a "mixed" expectation value between \&y and ^(r): 



(0(r)) A/ = 



(V(t)\Q\* v ) 
(*(r)|*y) : 



(6) 



where the operator O acts on the trial function vffy. The 
desired expectation values would, of course, have ^(r) on 
both sides; by writing vE'(r) = *y + <5\I>(t) and neglect- 
ing terms of order [(5^(r)] 2 , we obtain the approximate 
expression 



(0(r)> = 



(^(r)|0|^(r)) 
(*(r)|*(r)> 
(0(t)) m + [(0(r)) M 



(Oh 
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where (O)v is the variational expectation value. 



For the energy, the mixed estimate of Eq. (j6]) with O = 
W is itself a strict upper bound to the ground state for 
the simpler Hamiltonian, as can be seen by commuting 
half the imaginary time operator from the left to right 
hand side, giving 



(H'(t)) m = 



^{r/2)\H^{r/2)) 
(*(r/2)|*(r/2)) ' 



(8) 



The total energy is then given by this mixed estimate 
for H' plus the small difference ((H — H')) evaluated by 
Eq.©. 

For off-diagonal matrix elements required by transi- 
tions the generalized mixed estimate is given by the ex- 
pression 

(*f(T)\0\*'(T)) 



(0(r)) Mi + (0{t)) M} - (O)v , (9) 



where 
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and (0(T))nji is defined similarly. For more details 
see Eqs. (19-24) and the accompanying discussions in 

Rcf. g. 

Sources of systematic error in the GFMC evaluation 
of operator expectation values (other than H') include 
the use of mixed estimates and the constrained path al- 
gorithm for controlling the Fermion sign problem in the 
propagation of x P(t). These are discussed in Ref. [lj|; 
the convergence of the current calculations is addressed 
at the beginning of Sec. IIV1 



III. THE ELECTROMAGNETIC CURRENT 
OPERATOR 

The nuclear EM current operator is expressed as an 
expansion in many-body terms. The current utilized in 
this work includes up to two-body terms. In what follows, 
we use the notation 



k=(ki-k 2 )/2 , 



K i = (p<+ Pi )/2 , (12) 
K = K! + K 2 , (13) 



where pi (pj) is the initial (final) momentum of nucleon 
i, and q = ki + k 2 is the momentum associated with the 
external EM field. 

The one-body operator at leading order — or impulse 
approximation (IA) operator — consists of the convection 
and the spin-magnetization currents associated with an 
individual nucleon. It is derived from the non-relativistic 
reduction of the covariant single-nucleon current by ex- 
panding it in powers of pj/mjv, where tun is the nucleon 
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FIG. 1: Diagrams illustrating one- and two-body xEFT EM 
currents entering at LO (eQ~ 2 ), NLO (eQ _1 ), N2LO (eQ°), 
andN3LO (eQ 1 ). Nucleons, pions, and photons are denoted 
by solid, dashed, and wavy lines, respectively. 



mass and retaining the leading-order term. It reads 

[2ejv,iKi +inN,ivi x q] , (14) 



•IA 



2 ton 



where 



eN = {l + T z )/2, k n = (ks + kv t z )/2 , HN = e N + n N . 

(15) 

Here Kg = —0.12 n.m. and ny = 3.706 n.m. are the 
isoscalar (IS) and isovector (IV) combinations of the 
anomalous m.m.'s of the proton and neutron, e is the elec- 
tric charge, and t z is the Pauli isospin projection equal 
to +1 for protons and —1 for neutrons. 

The calculations of the m.m.'s of the A < 9 nuclei have 
been carried out utilizing two models for the two-body 
EM current operator, which are discussed in the next two 
subsections. 



A. xEFT current operator 

Two-body EM currents have been derived in re- 
cent yea rs with in p ionful chiral effective field formula- 
tions [lJMl20|-|22| (for a comparison of the different for- 
malisms we refer to the last four cited references). Here, 
we utilize the operators constructed within the formal- 
ism developed in Refs. [TlL fl^|. The xEFT operators are 
expanded in powers of pions' and nucleons' momenta, 
and consist of long- and intermediate-range components 
which are described in terms of one- and two-pion ex- 
change contributions, as well as contact currents encod- 
ing short-range mechanisms unresolved at the given or- 
der. These operators involve a number of low energy 
constants (LECs) which are fixed to the experimental 
data. The operators and fitting procedure have been re- 
cently described in Ref. [l3| . We refer to that work for a 
complete listing of the operators utilized in the present 



TABLE I: Dimensionless values of the isoscalar and isovector 
LECs entering the xEFT current operator at N3LO corre- 
sponding to cutoff A = 600 MeV. 



A 


C[ 5 x A 4 


4 x A 2 


C 16 x A 4 


d' s x A 2 


600 


5.238 


-0.2033 


-1.025 


4.980 



calculations, and limit ourselves to discussing the vari- 
ous contributions and to summarizing the fitting strategy 
adopted to constrain the LECs. 

The xEFT EM operators are diagrammatically repre- 
sented in Fig. [TJ They are expressed as an expansion 
in the low- momentum scale Q. Referring to Fig. [I] the 
leading-order (LO) term is counted as eQ -2 , and corre- 
sponds to the IA one-body operator given in Eq. (|14[) . 
The NLO term (of order eQ" 1 ) consists of the seag- 
ull and pion- in- flight one-pion-exchange (OPE) currents. 
These purely isovector currents involve two known LECs, 
the axial coupling constant gA — 1.29, and the pion de- 
cay amplitude = 184.6 MeV. The value for gA is 
determined from the Goldberger-Treiman relation gA = 
9-kNN Ft T /(2 itin), where the irNN coupling constant is 
taken to have the value gl NN /{A-K) = 13.63±0.20 [HHl]. 
The N2LO one-body contribution (of order eQ°) is a 
relativistic correction to the IA operator, and is thus ex- 
pressed in terms of the nucleons' experimental m.m.'s. 

At N3LO (e Q) we distinguish among four kinds of cur- 
rents. The first one ("LOOP" in the tables), accounts 
for the one-loop contributions of diagrams (e)-(i) and 
(l)-(o). These terms lead to a purely isovector current 
which involves the known LECs gA and . For diagram 
(m) we use the expression given in Ref. 1 1 31] which differs 
from that given in previous works [ll|j Il2l| by some of the 
present authors. 

Next we account for the contact currents illustrated in 
panel (k). We distinguish between minimal (MIN) and 
non-minimal (NM) currents. The former is linked to the 
vEFT contact potential at N2LO via current conserva- 
tion; therefore it involves the same LECs entering the 
NN potential, and is [U, El 
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where the low-energy constants Ci,...,CV, have been 
constrained by fitting np and pp elastic scattering data 
and the deuteron binding energy. We take their values 
from the Machlcidt and Entem 2011 review paper [25j |. 
Unknown EM LECs enter the NM current at N3LO 
which is given by 
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and the determination of the LECs C 15 and C 16 is dis- 
cussed below. 

The N3LO OPE current, diagram (j) in Fig. [H is given 
by Bill 
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We fix the LECs multiplying the isovector operators by 
relating them, in a resonance saturation picture, to the 
couplings of the TV to A excitation, i.e., 



/x 7 jva h>A 



9 mjv (m a — tun ) 



d 2i ~ 4 



(19) 



where toa and Tiyi are the mass of the A (hia — niN — 
294 MeV) and the TV to A axial coupling constant, 
M7JVA = 3 n.m. is the transition magnetic moment [26| . 
and h A /F^ = Una/iti^ with f% NA /(4rt) ~ 0.35 as ob- 
tained by equating the first-order expression of the A- 
decay width with the experimental value. The current 
proportional to these LECs, d' s and d 21 , reduces to the 
conventional TV- A current. The isoscalar term in Eq. (1181) 
saturates the standard pirj transition current fl2l. [l3l |20| . 

The EM operators described above have power-law be- 
havior for large momenta and need to be regularized be- 
fore they can be inserted between nuclear wave functions. 
The regularization procedure is implemented by means of 
a cutoff of the form 1131 



C A (k) = exp(-fc 4 /A 4 ) 



(20) 



We utilize the %EFT operators within a hybrid con- 
text, in which the matrix elements are evaluated with 
wave functions that are solutions of the realistic Hamil- 
tonian given in Eq. Intrinsic to this approach is a 
mismatch between the short-range behavior of the nu- 
clear potential and that of the EM operator. As a con- 
sequence, the current is not strictly conserved. 

The fitting of the unknown LECs entering the EM cur- 
rents, namely C' 15 , C' w , d' s , and d' g (with d' 21 = dg/4, 
as implied by the A-saturation mechanism) was done in 
Ref. [1J|. In that work, the cut-off A was varied in the 
range (500-600) MeV and the LECs were constrained 
to reproduce a set of nuclear EM observables for any 
given A in this range. Three different parametrizations 
were tested in the A = 2-3 nuclei. The trinuclcon wave 
functions, required for the evaluation of the matrix ele- 
ments, were obtained with the hyperspherical harmonics 
(HH) expansion discussed in Refs. |27H29| with a nuclear 
Hamiltonian consisting of the Argonne V\% (AV18) i and 
Urbana IX (UIX) [U potentials. (There is very little dif- 
ference between the A = 3 wave functions for AV18+UIX 
and AV18+IL7.) The three models (labeled model I, II, 
and III) determine the LECs multiplying isoscalar oper- 
ators (i.e. C 15 and d' 9 ) so as to reproduce the experi- 
mental dcutcron m.m. and the isoscalar combination of 



the trinucleon m.m.'s. Models II and III fix the isovector 
LEC d' s by A-saturation as indicated in Eq. (fT9")l . while 
C 16 is constrained so as to reproduce either the np radia- 
tive capture cross section at thermal neutron energies in 
model II, or the isovector combination of the trinucleon 
magnetic moments in model III. In model I, d' s is left as a 
free parameter and is constrained, along with C{ 6 , so as 
to reproduce both the np radiative capture cross section 
and the isovector combination of the trinucleon m.m.'s. 
As already observed in Ref. [li|, model I leads to un- 
naturally large values for both isovector LECs, severely 
spoiling the convergence pattern of the chiral expansion. 
We have, nevertheless, tested all three models (with cut- 
offs of both 500 and 600 MeV) in VMC m.m. calculations 
for A = 3-8 nuclei, and verified that this pathology — i.e., 
the lack of convergence — persists, and indeed gets worse, 
in larger systems. We have therefore disregarded model 
I, and adopted model III with cut-off A = 600 MeV in 
the present study. The parameters entering this model, 
obtained from the calculations performed in Ref. [l3j , are 
listed in Table U Model III (with A = 600 MeV), when 
tested in VMC calculations, produced the best results for 
the m.m.'s. 



B. SNPA current operator 

The two-body currents in the SNPA formalism have 
been described in detail most recently in Ref. [l(| ■ These 
currents are separated into model-independent (MI) and 
model-dependent (MD) terms. The former (MI) are de- 
rived from the TVTV potential (the AV18 in present case), 
and their longitudinal components satisfy, by construc- 
tion, current conservation with it. They contain no free 
parameters, and their short-range behavior is consistent 
with that of the potential. The dominant terms, isovector 
in character, originate from the static part of the poten- 
tial, which is assumed to be due to exchanges of effective 
pscudoscalar (PS or "7r-like") and vector (V or "p-likc") 
mesons. The associated currents are then constructed 
by using the PS and V propagators, projected out of 
the static potential [l(|. Additional (short-range) cur- 
rents follow by minimal substitution in the momentum- 
dependent part of the potential. They have both isoscalar 
and isovector terms, and lead to contributions which are 
typically much smaller (in magnitude) than those gener- 
ated by the PS and V currents. At large inter-nucleon 
separations, where the TVTV potential is driven by the 
OPE mechanism, the MI current coincides with the stan- 
dard seagull and pion-in-flight OPE currents diagram- 
matically illustrated in panels (b) and (c), respectively, 

of Fig. n 

The MD currents are purely transverse, and uncon- 
strained by current conservation. The dominant term 
is associated with excitation of intermediate A isobars, 
which arc treated non-perturbativcly with the transition 
correlation-operator method developed in Ref. 31 



These A currents arc discussed in considerable detail in 



5 



Ref. 0] ■ Additional (and numerically small) MD currents 
arise from the isoscalar pirj and isovector ujtt"/ transition 
mechanisms. The values for the coupling constants en- 
tering them are also listed in Ref. Q. 



IV. RESULTS 

The IA m.m. for the 8- and 9-body nuclei have signif- 
icantly higher Monte Carlo statistical errors than most 
quantities that we have computed with GFMC. There- 
fore we present two examples of the GFMC propagation 
as a function of the imaginary time r. Fig. [2] shows the 
propagation of a typical case, 8 Li. Three propagations 
are shown, one in which the constrained propagation [l9j 
is relaxed with 40 unconstrained steps (n u = 40) and two 
with n u = 80. The n u = 40 case was made with 10,000 
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FIG. 2: (Color online) Propagation of the energy (a) and 
m.m. (b) as a function of imaginary time r for the ground 
state of 8 Li. 



the starting VMC values are within 5% of the final re- 
sults. As has been shown in Ref. these quantities 
are all converged by n u = 40, often by n u = 20. The 
quadrupole moments are much more difficult to evaluate, 
because they have long-lived oscillations in the propaga- 
tion time r. 

The IA and %EFT MEC m.m. are shown in panel (b). 
The statistical fluctuations of the IA term are much worse 
than those of the MEC term when the same number of 
configurations are used. Also there may be a small sys- 
tematic change in the IA term going from n u = 40 to 
n u = 80; the average values are 1.120(27) and 1.164(17), 
respectively, giving a difference of 0.044(32). The MEC 
does not have this sensitivity; the two calculations are 
in excellent agreement. However because of the rapid 
growth of statistical error with increasing n u , the n u = 80 
calculation needs seven times as many walkers to achieve 
the same statistical error. 

The Monte Carlo statistical errors in our computed 
m.m.'s are most severe for 9 C, as is shown in Fig. [3] Two 
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walkers while the n u = 80 cases are averages of calcula- 
tions with 20,000 and 40,000 (60,000 total) and 20,000, 
40,000 and 80,000 walkers (140,000 total), respectively. 
The energy is shown in panel (a); it is similar to results 
shown in Ref. [l9| . As can be seen, there is a rapid drop 
from the initial VMC value at r = that reaches a sta- 
ble result before 0.1 MeV . The results for all quantities 
presented in this article are averages over r from 0.2 to 
0.8 MeV -1 , as indicated by the solid lines, with statisti- 
cal errors denoted by the dashed lines. The propagations 
of point proton and neutron radii are similar, except that 



FIG. 3: (Color online) Propagation of the energy (a) and 
m.m. (b) as a function of imaginary time r for the ground 
state of 9 C. 

propagations are shown, one in with n u = 60 and one 
with n u = 80. The first was made using 20,000 walkers 
while the latter is the average of two calculations, each 
with 160,000 walkers. Again the energy, shown in the 
panel (a), is well converged for both cases. As before for 
8 Li, the xEFT MEC m.m., panel (b), has smaller fluctu- 
ations for a given number of walkers than the IA m.m. 
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In this case the xEFT MEC is a large (80%) addition to 
the IA m.m. 

In both examples the statistical fluctuations in the IA 
term are much larger than in the MEC term. However 
the evaluation of the MEC requires much more compu- 
tational effort per walker than does the propagation and 
IA term. Therefore for most of the calculations, we prop- 
agate a large number of walkers using n u = 60 or 80 to 
obtain the IA m.m. term (and also the other reported 
quantities). The MEC is obtained with comparable sta- 
tistical error using fewer walkers and the two numbers 
and their errors combined to get the total m.m. The 
propagations are averaged over r = 0.2 to 0.8 McV -1 . 

The large statistical fluctuations (and possible n u sen- 
sitivity) are coming from the IV combination of the spin 
term of the IA m.m., Eq. (|T4|) . The IV convection term 
and both IS terms have much smaller fluctuations. Thus 
if isospin symmetry is assumed for the wave functions of 
isobaric analogs, we can make precise statements about 
the IS m.m. However if we do not want to assume such 
isospin symmetry, as in the 9 C- 9 Li case below, then we 
have to make separate calculations for each nucleus and 
the large errors in the IV parts make the extraction of an 
IS m.m. with small statistical error impossible. 

The energies E, point proton rms radii r p (and point 
neutron rms radii r n for N ^ Z nuclei), m.m.'s fj, in IA, 
and quadrupole moments Q for the nuclear states calcu- 
lated in this work are presented in Table HIl along with ex- 
perimental values where available. Experimental energies 
arc from Ref. [33|, EM moments are from Refs. [34j - l37j ]. 
and point radii are converted from the charge radii given 
in Refs. [38M42l |. Many energies for A < 7 nuclei evalu- 
ated with the AV18+IL7 Hamiltonian have been reported 
previously in Ref. f43j . The present energies, which are 
from independent calculations, are in agreement with the 
previous results within the Monte Carlo statistical errors 
shown in parentheses. 

For many of the isobaric analog states, the energy 
and moments are calculated using the GFMC wave func- 
tions generated for the T z = — T state and then sim- 
ply interchanging protons and neutrons to evaluate the 
T z = +T state. These calculations are denoted by an as- 
terisk (*) in the table, and will be referred to as charge- 
symmetry-conserving (CSC) results. For 3 He, 8 B, and 
9 C ground states we also made independent calculations 
with different starting VMC wave functions and different 
isoscalar Coulomb terms [J, |32| in the GFMC propaga- 
tor appropriate to the T z = +T state. We then use these 
wave functions to predict the quantities in their isobaric 
analogs, i.e., 3 H, 8 Li, and 9 Li. The pairs of independent 
solutions for the isobaric analogs will be referred to as 
charge-symmetry-breaking (CSB) results. Thus six nu- 
clear states have two entries in the table, comparing a 
direct calculation with the prediction by charge symme- 
try from its isobaric analog. 

We see from the table that the energies in these paired 
calculations are in generally good agreement, with the 
largest discrepancy for A = 9, where the difference is 



~ 2% and the statistical errors almost touch. In each of 
the CSB cases, the T = — T z state is more bound than its 
isobaric analog, and the expectation values of individual 
terms in the nuclear Hamiltonian, like {Ki) and are 
larger in magnitude. 

The point nucleon rms radii are slightly larger for the 
proton-rich nuclei compared to the charge symmetric so- 
lution from the proton-poor isobaric analog. For exam- 
ple, the proton rms radius of 9 Li ( 9 C) is smaller (larger) 
when the appropriate Coulomb term is included in the 
GFMC propagator, indicating that the system is more 
compact (diffuse). If 9 Li ( 9 C) is constructed from the 
9 C ( 9 Li) solution, then it appears to be a more diffuse 
(compact) system. Consistently with this weak sensitiv- 
ity of the calculated energies and radii to the charge sym- 
metry picture implemented to derive the nuclear wave 
functions, we find that the calculated m.m.'s in IA are 
not statistically different in the T = h , 1 cases and we 
see only very weak evidence that the IA m.m.'s of the 
A = 9, T = 3/2 systems are affected by charge symme- 
try. The quadrupole moments are also fairly consistent 
in the paired results and close to the experimental values. 

A. Magnetic Moments in A=2— 9 Nuclei 

The calculations of the matrix elements, both diagonal 
and off diagonal, have been described in detail in Refs. 0, 
0|. In particular, the I A matrix element is evaluated 
using the Ml operator induced by the one-body current 
given in Eq. (|14[) . namely 

M IA = ^2(e N ,iLi +hn,i(Ti) ■ (21) 

i 

The matrix element associated with the MEC contribu- 
tion is 

(JJ,M f | M MEC | Jf, Mi) = 

- nim ^L {r M \jf^{ q ±)\ jr, Mi ) , (22) 

q— S-0 ? 

where the spin-quantization axis and momentum transfer 
q are, respectively, along the z and x axes, and Mj = J. 
The various contributions are evaluated for two small val- 
ues of q < 0.02 fm _1 and then extrapolated smoothly to 
the limit q=0. The error due to extrapolation is much 
smaller than the statistical error in the Monte Carlo sam- 
pling. 



In Table UTTl we show, in addition to the proton and neu- 
tron experimental m.m.'s, the experimental and calcu- 
lated m.m.'s for the A = 2 and 3 nuclei, including MEC 
contributions from the EM currents in the SNPA and 
%EFT models. In the table we label with IS and IV the 
isoscalar and isovector combinations of the magnetic mo- 
ments as given by: 

fi{T,T z ) = fi(TS) + fi(JV)T a . (23) 
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TABLE II: GFMC results for A < 9 nuclear states studied in this work, compared to experimental values [33H43] . Numbers in 



parentheses are 


statistical errors 


for the GFMC calculations or 


experimental 


errors; errors 


of less than 


one in the last decimal 


place are not shown. 
















A Z(J*,T) 


E [MeV 




r p [r n ] [fm] 




//(IA) [n.m.] 


Q [far 






GFMC 


IZjXpt. 




Expt. 


GFMC 


Expt. 




T7 , . ... 1 
JLXpt . 


X H(1+, 0) 


-2.225 


—2.2246 


1.968 


i.y /d(oj 


0.847 


0.8574 


0.270 


0.286 


3tjv 1 + 1 \ 
H (2 .5) 


—8.50(1) 


—8.482 


1 58h 761 


1 c; Q ^1 n\ 
1 .OO^IUJ 


2.556 


2.979 






5 TI / 1 -4- 1 \ J- 

3 H(§ ,|)* 


-8.46(1) 




1.60[1.80] 




2.550 








3He(| + ,I) 


-7.73(1) 


—7.718 


1 80h 601 


1 Vfi/I \ 

t-tb(L) 


-1.743 


-2.127 






3 Hc(i + ,i)* 


-7.75(1) 




1.76[1.58] 




— 1. / OU 








6 Li(l+,0) 


-31.82(3) 


—31.99 


2.39 


2.45(4) 


0.817 


0.822 


—0.20(6) 


—0.082(2) 


6 Li(0+,l) 


-28.44(4) 


—28.43 














'-"■\ 2 > 2' 


—39.0(1) 


-39.24 


2.28[2.47] 


2.31(5) 


2.87 


3.256 


-4.0(1) 


-4.00(3) 




—38.9(1) 


-38.76 














7 Bc(| ,±)* 


-37.4(1) 


-37.60 


2.47[2.28] 


2.51(2) 


-1.06 


-1.398(15) 


—6.7(1) 




7 Bc(i 


-37.3(1) 


—37.17 














8t \( r >+ 1 ^ 




-41.28 


2.10[2.46] 


2 20f5"l 


1.16(2) 


1.654 


3.3(1) 


3.14(2) 


8 Li(2+, 1)* 


-41.0(2) 




2.11[2.48] 




L.Loyo ) 




3.0(4) 




8 Li(l+,l) 


-40.1(2) 


-40.30 














8 Li(3+,l) 


-38.5(3) 


-39.02 














8 B(2+,1) 
°B(2^, 1)* 


-37.5(2) 


—37.74 


2.48[2.11J 




1.45(1) 


1.036 


5.9(4) 


n 01/01 \ 

6.83(21) 


—37.8(2) 




2.46[2.10] 




1.42(2) 




6.5(2) 




8 BC1 + 1"!* 


—36.6(2) 


-36.97 














8 B(3+.l)* 


-34.8(3) 


-35.42 














9 Be(|-,i) 


-58.1(2) 


-58.16 


2.37(1)[2.56(1)] 


2.38(1) 


-1.18(1) 


-1.178 


5.1(1) 


5.29(4) 


9 Be(|-,I) 


-55.7(2) 


-55.74 














9 B(|-,|)* 


-56.3(1) 


-56.31 


2.56(1)[2.37(1)] 




2.97(1) 




4.0(3) 




9 B(|-,ir 


-53.9(2) 


-53.95 














9 Li(|M) 


-45.2(3) 


-45.34 


1.97(1)[2.36(2)] 


2.11(5) 


2.66(3) 


3.437 


-2.3(1) 


-3.06(2) 


9 Li(§ ,§)* 


-45.9(3) 




2.03(1)[2.45(1)] 




2.64(4) 




-2.7(2) 






-43.2(4) 


-42.65 














9 c(|-,|) 


-39.7(3) 


-39.04 


2.45(1)[2.03(1)] 




-0.75(3) 


-1.391 


-4.1(4) 




9 c(|-,|r 


-38.8(3) 




2.36(2)[1.97(1)] 




-0.82(4) 




-3.7(1) 





With the label MEC we denote anything that goes be- 
yond the IA picture, therefore the xEFT MEC current 
includes also the one-body relativistic correction opera- 
tor entering at N2LO. The results for the deuteron are 
from calculations of matrix elements with wave functions 
which are exact solutions of the two-body Schrodingcr 
equation with the AV18 potential. Results for the A = 3 
nuclei are from GFMC calculations with the nuclear 
Hamiltonian consisting of the AV18 two-body and IL7 
three-body potentials (AV18+IL7), while those desig- 
nated with HH are results from hyperspherical harmon- 
ics from Ref. [l3| obtained with the nuclear Hamilto- 
nian consisting of the AV18 and the UIX three-body po- 
tentials (AV18+UIX). Both the GFMC and HH wave 
functions have been constructed separately without ex- 
ploiting charge symmetry. Strictly speaking, the GFMC 
m.m.'s are for the propagating Hamiltonian H', i.e., 
AV8'+IL7', as discussed in Sec. HI The small 0.3%-0.5% 
difference between the GFMC and HH IA values may be 
attributable to this difference in the Hamiltonians. 

Numerical differences between the calculated %EFT 
MEC terms are also affected by an additional approx- 
imation implemented in the GFMC calculations. The 
one-body m.m. operator associated with the relativistic 



correction at N2LO (illustrated in panel (d) of Fig. [T]) 
reads lllll 



N2LO 



E 



"N , := i 



{pi , ejy.iLj + i-i N ,iCn} 



-eN,iPi x (<Ti x 



(24) 



where Pi = —iVi and L, are the linear momentum and 
angular momentum operators of particle i, and {. . . , . . .} 
denotes the anticommutator. In the GFMC calculations 
we do not explicitly evaluate the pf term, but instead ap- 
proximate it with its average value, that is pf ^< pf >, 
as determined from the expectation value of the kinetic 
energy operator in each nucleus. This approximation 
leads, in the case of 3 He, to a 5% difference in the MEC 
correction which itself is a 20% correction to the total 
calculated m.m., and to an even smaller effect in 3 H. 
The values utilized in the GFMC calculations are re- 
ported in Table IIV1 for the nuclei investigated in this 
work. Of course, this approximation only affects the 
N2LO xEFT operator. The HH calculations use the non- 
approximated operator at N2LO and are designated with 
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TABLE III: Magnetic moments in nuclear magnetons for A = 2-3 nuclei evaluated with SNPA and xEFT EM current operators. 
The current model labeled xEFT^ accounts for an exact calculation of the N2L0 relativistic correction to the IA current (see text 
for explanation). Results labeled with HH are obtained utilizing trinucleon wave functions constructed with the hyperspherical 
harmonics (HH) method developed in Refs. [27l - [2^ |. and a nuclear Hamiltonian consisting of AV18+UIX. The remaining results 
(except for the deuteron ones which are exact) are from GFMC calculations discussed in the text. The spatial symmetry (s.s.) 
of the nuclear wave function is also given. The IS and IV labels indicate the isoscalar and isovector combinations. Results 
obtained with the ^EFT and xEFT"^ models are not predictions (see text for explanation). 



Nucleus (J*";T) Current s.s. IA MEC Total Expt. 



n(I + ;i) 

Kl + ;§) 
is 

IV 



-1.913 
2.793 
0.440 
4.706 



2 H(| + ;i) SNPA [2] 0.8470 

XEFT 0.8470 
xEFT f 0.8472 



0.0012 0.8482 0.8574 

0.0134 0.8604(1) 
0.0102 0.8574 



SNPA 




[3] 


2.556(1) 


XEFT 






2.556(1) 


xEFTt 


(HH) 




2.569 


SNPA 




[3] 


-1.743(1) 


XEFT 






-1.743(1) 


xEFT f 


(HH) 




-1.749 


SNPA 






0.407 


X EFT 






0.407 


xEFT f 


(HH) 




0.410 


SNPA 






-4.299 


XEFT 






-4.299 


XEFT 1 


(HH) 




-4.318 



0.347(2) 


2.903(2) 


2.979 


0.404(1) 


2.960(1) 




0.410 


2.979 




-0.334(2) 


-2.077(2) 


-2.127 


-0.357(1) 


-2.100(1) 




-0.378 


-2.127 




0.006 


0.413 


0.426 


0.024 


0.431 




0.016 


0.426 




-0.681 


-4.980 


-5.106 


-0.761 


-5.060 




-0.788 


-5.106 





TABLE IV: Values (in units of fm 2 ) of the approximated 
{pi), entering the xEFT current at N2LO, used in the GFMC 
calculations. 



A 


2 3 


6 7 8 9 


(pi) 


0.5 0.8 


1.2 1.2 1.2 1.3 



XEFTt. These calculat ions have been used to constrain 
the LECs entering the xEFT currents, therefore the HH 
results presented in Table Hill arc not predictions, in that 
they reproduce the experimental data by construction. 



There are also tiny numerical differences between the 
SNPA calculations presented here and those reported in 
Ref. These may be due to differences in the starting 
variational wave functions as well as systematic uncer- 
tainties in the GFMC calculation. At any rate, all the 
numerical differences mentioned above are of little impor- 
tance if one accounts for the sensitivity of these results to 
the nuclear and current models utilized. The results in 
the A = 2 and 3 nuclei show that the SNPA model under- 
estimates the isoscalar component in both the deuteron 



and the trinucleon m.m.'s. 



In Table [V] we report the GFMC calculations for the 
m.m.'s of the A = 6-9 nuclei. We compare results ob- 
tained using either the SNPA or the xEFT MEC cur- 
rents. The MEC corrections evaluated in both models 
are qualitatively in agreement. They boost the IA in the 
direction of the experimental data in all cases, except 
for 6 Li and 9 Be. In these systems the IA results are al- 
ready in very good agreement with the experimental data 
while the small MEC contributions make the predictions 
slightly worse. 

The calculations denoted with an asterisk (*) are ob- 
tained exploiting charge symmetry, that is by interchang- 
ing protons and neutrons to generate the isobaric analogs 
wave functions. For the A = 7,8 nuclei, we verified 
that m.m. predictions obtained with independent nuclear 
wave functions for the isobaric analogs, i.e., calculations 
which account for the proper isoscalar Coulomb term 
in the starting VMC wave functions as well as in the 
GFMC propagator, are essentially identical to the stan- 
dard charge symmetric results, and therefore we do not 
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TABLE V: Magnetic moments in nuclear magnetons for A = 6-9 nuclei evaluated with SNPA and xEFT EM current operators. 
Results labeled with a star are obtained exploiting charge symmetry. The dominant spatial symmetry (s.s.) of the nuclear 
wave function is given. 



Nucleus(J' r ;r) Current 


s.s. 


IA 


MEC 


Total 


Expt. 




SNPA 


[421 


0.817(1) 


-0.010(1) 


0.807(1) 


0.822 








0.817(1) 


0.020(1) 


0.837(1) 




7t • / 3 — 1 \ 

Li(| ;§) 


SNPA 


[43J 


O Q7/1 \ 
Z.O 


V.Zo{Z ) 


o.L2{2 ) 


O.ZOD 




XEl 1 




ij.Ol 1 -L I 


U.O i y J- j 


^ 9d( 1 ^1 

O.^j'I^ J. ^ 




7 Be(r;i)* 


SNPA 


[43] 


-1.06(1) 


-0.39(2) 


-1.45(2) 


-1.398 




X EFT 




— 1.06(1 J 


— 0.36(1) 


-1.42(1) 




IS 


SNPA 




0.90 


-0.07 


0.83 


0.929 




X EFT 




0.90 


0.01 


0.91 




IV 


SNPA 




-3.93 


-0.64 


-4.57 


-4.654 




X EFT 




—3.93 


—0.73 


—4.66 




Li(2 + ;1) 


SNPA 


I A Q 1 1 

[431] 


1.10(Z J 


U.2U(2 J 


l.O0(O J 


i a k a 
1.004 




V F,FT 




1.16(2) 


0.33(1) 


1.49(2) 




B(2 ,1) 


olN rA 


|4dlJ 


1.42(2) 


— 0.42(2) 


1.00(3) 


1.036 




X^b 1 




1.42(2) 


-0.31(1) 


1 11(2) 




IS 


SNPA 




1 9Q 


nil 

— U. 1 1 


IIS 
l.lo 


1.04:0 




XEFT 




1.29 


n m 


1 ^0 
i .ou 




IV 


SNPA 




H 1 ^ 
U. lo 


— u.oi 


— u. lo 


— u.ouy 




X-fji 1 1 




0.13 


—0.32 


—0.19 






QMP4 


[4^191 
[4JZJ 








0.4:0 J 




vEFT 




2.66(3) 


0.70(2) 


3.36(4) 






olNrA 


f4^9l 
[4JZJ 






1 9Q/'c;~\ 


— i.oyi 




Xhib 1 




-0.75(3) 


-0.60(3) 


-1.35(4) 




is 


SNPA 




u.yo 


— U.U ( 


u.oy 


i.uzo 




X EFT 




u.yu 


o 


i m 




T T r 

IV 


C\TT) \ 




11/1 
— 1.14 


n 97 


111 
— 1.41 


— i.ouy 




XEFT 




— 1.14 


—0.43 


— 1.57 




9 Be(| + ;i) 


SNPA 


[441] 


-1.18(1) 


-0.12(1) 


-1.30(1) 


-1.178 


X EFT 




-1.18(1) 


-0.11(1) 


-1.29(1) 




9 B(| + ;I)* 


SNPA 


[441] 


2.97(1) 


-0.10(1) 


2.87(1) 


n.a. 


X EFT 




2.97(1) 


0.09(1) 


3.06(1) 




IS 


SNPA 




0.89 


-0.11 


0.78 






XEFT 




0.89 


-0.01 


0.88 




IV 


SNPA 




4.15 


0.02 


4.17 






XEFT 




4.15 


0.20 


4.35 





report them. In the A = 9, T = | case there is weak evi- 
dence for a CSB effect so in this case we show the results 
of the two independent calculations. 

Also in the A — 6-9 nuclear m.m.'s, the difference 
between the SNPA and xEFT corrections is more pro- 
nounced in the isoscalar component. In all cases, the 
XEFT corrections are more positive (or less negative) 
than the corresponding SNPA. This makes the xEFT 
predictions closer to the experimental values. The isovec- 
tor corrections evaluated with the two models are rea- 
sonably in agreement with each other, although they are 
bigger when derived from the xEFT model. MEC correc- 
tions are crucial to bring the theory closer to the exper- 
imental values. Their effect is particularly pronounced 
in the isovector combination of the A = 9, T = 3/2 
nuclei's m.m.'s, for which the MEC SNPA (xEFT) cor- 
rection provides ~ 20% (~ 30%) of the total calculated 



isovector contribution. 

It is interesting to note that, despite the large effect 
observed in the A = 9, T = 3/2 systems, MEC correc- 
tions are considerably smaller in the ^4 = 9,T=l/2 
nuclei. This feature can be explained by considering the 
dominant spatial symmetry (s.s.) of the wave functions 
associated with the A = 9 systems. In particular, the 
dominant spatial symmetry of 9 Be ( 9 B) is [441], corre- 
sponding to an [a, a, n(p)] structure as shown in Ref. [4~4| . 
A single nucleon outside an a particle feels no net OPE 
potential, and this holds true also for a single nucleon 
outside a double-a [44] symmetry state. Consequently, 
the NLO OPE currents illustrated in panels (b) and (c) 
of Fig. [TJ which are generally the largest MEC terms 
in both SNPA and xEFT approaches, do not contribute 
significantly. On the other hand, the dominant spatial 
symmetry of 9 C ( 9 Li) is [432] ~ [a, 3 He ( 3 H), pp (nn)], 
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FIG. 4: (Color online) Magnetic moments in nuclear magne- 
tons for A < 9 nuclei. Black stars indicate the experimen- 
tal values |!35r:37j], while blue dots (red diamonds) represent 
GFMC calculations which include the IA one-body EM cur- 
rent (total xEFT current up to N3LO). Predictions are for 
nuclei with A > 3. 



and the NLO OPE term contributes in both the trinu- 
cleon clusters and in between the trinucleon clusters and 
the valence pp (nn) pair. The IA m.m. for 9 Bc is close 
to the experimental value, while those for 9 Li and 9 C 
are far from the data, so this pattern of small and large 
MEC corrections provides good overall agreement with 
the data. 

The xEFT results reported in Tables HIT] and IVl ar e 
summarized in Fig. 2J where the experimental data [3J- 
H3| (there are no data for the m.m. of 9 B) are repre- 
sented by black stars. We show also the experimen- 
tal values for the proton and neutron m.m.'s, as well 
as their sum, which corresponds to the m.m. of an S- 
wave deuteron. The experimental values of the A = 2-3 
m.m.'s have been utilized to fix the LECs, therefore pre- 
dictions are for A > 3 nuclei. The blue dots labeled 
as GFMC(IA) represent theoretical predictions obtained 
with the standard IA one-nucleon EM current entering 
at LO: diagram (a) of Fig. [TJ The GFMC(IA) results 
reproduce the bulk properties of the m.m.'s of the light 
nuclei considered here. In particular, we can recognize 
three classes of nuclei with non-zero m.m.'s, i.e., odd- 
even nuclei whose m.m.'s are driven by an unpaired va- 
lence proton, even-odd nuclei driven by an unpaired va- 
lence neutron, and odd-odd nuclei with either a deuteron 
cluster or a triton-neutron ( 3 Hc-proton) cluster outside 
an even-even core. Predictions which include all the con- 
tributions to the N3LO xEFT EM currents illustrated 
in Fig. [JJ are represented by the red diamonds of Fig. HI 
labeled GFMC(TOT). In all cases except 6 Li and 9 Be 
(where the IA is already very good and the MEC correc- 
tion is very small) the predicted m.m.'s are closer to the 
experimental data when the MEC corrections arc added 
to the IA one-body EM operator. 



It is also interesting to consider the spatial distribution 
of the various contributions to the m.m., i.e., to examine 
the magnetic density. The one-body IA contributions 
from the starting VMC wave functions are shown in Fig. [5] 
for the isobaric analog pairs 7 Li- 7 Be, 8 Li- 8 B, and 9 Li- 
9 C. (The VMC values for the IA m.m.'s arc within a few 
% of the final GFMC values, so we expect their spatial 
distribution to be reasonably accurate.) In the figure, the 
red upward-pointing triangles are the contribution from 
the proton spin, p P [p P f(r) — p P i(r)], and similarly the blue 
downward-pointing triangles are the contribution from 
the neutron spin. The green diamonds are the proton 
orbital (convection current) contribution, and the black 
circles are the sum. The integrals of the black curves over 
d 3 r give the total m.m.'s of the nuclei in IA. 

For the neutron-rich lithium isotopes, there is one un- 
paired proton (embedded in a p-shell triton cluster) with 
essentially the same large positive contribution in all 
three cases. The proton orbital term is also everywhere 
positive, but relatively small. For 7 Li and 9 Li, the neu- 
trons are paired up, and give only a small contribution, 
so the total m.m. is close to the sum of the proton spin 
and orbital parts. However 8 Li has one unpaired neu- 
tron which acts against the proton and significantly re- 
duces the overall m.m. values. For the proton-rich iso- 
baric analogs, there is one unpaired neutron (embedded 
in a p-shcll 3 He cluster) with the same sizable negative 
contribution in all three cases. In 7 Be and 9 C, the pro- 
tons are paired up and give little net contribution, but 
the orbital term is always positive and acts against the 
neutron spin term. In 8 B there is also one unpaired pro- 
ton, which gives a bigger contribution than the unpaired 
neutron and results in a net positive m.m. value. 

In Table I VII we explicitly show the various contribu- 
tions entering the xEFT operator. The labeling in the 
table has been defined in Sec. IIII A[ We list the contribu- 
tions at each order. At N3LO, we separate the terms that 
do not depend on EM LECs (i.e. the LOOP contribution 
and the contact MIN currents; the former depends on the 
known axial coupling constant, gA, and pion decay am- 
plitude, Fa, while the latter depends on the strong LECs 
entering the AW xEFT potential at N2LO) and those 
that depend on them (i.e. the contact NM and the OPE 
current whose isovector component has been saturated 
with the A transition current). In most cases, chiral 
convergence is observed but for the isovector N3LO OPE 
contribution whose order of magnitude is in some cases 
comparable to the OPE contribution at NLO. It is likely 
that the explicit inclusion of A degrees of freedom in the 
present %EFT would significantly improve the conver- 
gence pattern, since in such a theory this isovector OPE 
current, presently entering at N3LO, would be promoted 
to N2LO. 

In Table I VII we do not provide the errors associated 
with the individual terms at each order because they are 
highly correlated. We limit ourselves to report the errors 
associated with the IA, MEC, and total results. Also 
in this table, we denote calculations performed enforcing 
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FIG. 5: (Color online) Magnetic density in nuclear magnetons per fm 3 for selected nuclei (see text for explanation). 



charge symmetry with an asterisk (*). In these calcula- 
tions, the isoscalar component of pure isovector operators 
(that is the OPE operator at NLO and the LOOP oper- 
ator at N3LO) are obviously zero. Calculations in which 
the nuclear wave functions are constructed independently 
present an isoscalar admixture even in the purely isovec- 
tor corrections, namely the NLO and N3LO LOOP. We 
do not report the individual contributions entering the 
SNPA calculations. For the A < 7 nuclei, they are found 
to be in agreement with those reported in Ref. 0- 

The A — 9, T = 3/2 nuclei are very interesting systems 
not only for the large effect produced by the isovector 
MEC correction, but also for the 'anomaly' associated 
with their isoscalar component. In Ref. [45| and refer- 
ences therein, the role of mirror symmetry — or charge 
symmetry — in the A = 9, T = 3/2 nuclei has been in- 
vestigated. In particular, if mirror symmetry is assumed, 
that is if the equality of nn and pp forces is enforced, 
then the isoscalar spin expectation value < a z > can be 
deduced as 



< cr z >= 



H(T Z = +T) + n(T z = -T)-J _ gggg) - J 

(5? + 5?-l)/2 ~ 0.3796 ' 

(25) 

In the equation above, /i and /^(IS) are the total m.m. 
and the isoscalar combination of the mirror nuclei m.m.'s, 

J is the total angular momentum and g p s ^ is the spin g 
factor of the proton (neutron). The experimental value 
of < a z > is 1.44 [46| which is — quoting from Ref. [45l]— 
"anomalously large if compared to the single-particle es- 



timate of 1." 

Driven by the discussion reported in Ref. (45| . we cal- 
culate the < a z > value for the A = 9,T = 3/2 nu- 
clei using the CSB and CSC models described at the 
beginning of this section. To make this last comment 
explicit, we show three different calculations for the 9 Li 
and 9 C m.m.'s in Table EU using the xEFT MEC. The 
first calculation is performed by constructing the nuclear 
wave functions independently including the appropriate 
Coulomb term in the GFMC propagator for each nucleus. 
The second (third) CSC calculation uses a 9 C ( 9 Li) wave 
function which is constructed from that of 9 Li ( 9 C) by 
interchanging protons with neutrons, i.e., by imposing 
mirror symmetry. The CSB isoscalar component is ob- 
tained combining the m.m.'s of 9 C and 9 Li obtained in- 
dependently. Therefore the error associated with this 
observable follows from propagating the statistical errors 
of the calculated m.m.'s of the A = 9, T = 3/2 nuclei. 
Within the CSC picture, the error associated with the 
isoscalar combination of the A = 9, T = 3/2 m.m.'s can 
be directly evaluated in the GFMC calculation. These 
CSC isoscalar values are found to be very stable, with a 
statistical error of less that ~ 1%. 

We use the CSC and CSB IA and total values for the 
isoscalar combination of the 9 C and 9 Li m.m.'s to eval- 
uate the < a z > value as given in Eq. (f2"5|) . The results 
are shown in Table IVIIII From this Table, we see that 
the CSC values arc consistent with the single particle 
prediction of 1. We observe that, in the CSB calculation 
both the isoscalar IA and MEC corrections are larger 



12 



TABLE VI: Magnetic moments in nuclear magnetons of A < 9 nuclei evaluated with the xEFT current operator. Individual 
contributions entering at LO (IA), NLO, N2LO, N3LO are shown (see text for explanation). Remaining notation is as defined 
in Table [V] 



Nucleus(J' r ;T) s.s. IA NLO N2LO N3LO N3LO-LECs MEC Total Expt. 



LOOP MIN NM OPE 

2 H(l+;0) [2] 0.8470 -0.0042 - 0.0364 -0.0135 -0.0052 0.0135 0.8604 0.8574 

3 H(i + ;i) [3] 2.556 (1) 0.253 -0.033 0.058 0.035 -0.011 0.102 0.404 2.960 (1) 2.979 

3 He(| + ;i) [3] -1.743 (1) -0.248 0.024 -0.055 0.056 -0.022 -0.110 -0.357 -2.100 (1) -2.127 

IS 0.407 (1) 0.002 -0.005 0.001 0.046 -0.017 -0.004 0.023 0.430 (1) 0.426 

IV 4.299 (1) -0.501 0.057 -0.113 0.020 -0.011 -0.213 -0.760 -5.060 (1) -5.106 

6 Li(l+;0) [42] 0.817 (1) - -0.012 0.063 -0.023 -0.007 0.020 (1) 0.837 (1) 0.822 

7 Li(|";i) [43] 2.87 (1) 0.237 -0.062 0.064 0.034 -0.012 0.107 0.37 (1) 3.24 (1) 3.256 

7 Be(|";i)* [43] -1.06 (1) -0.237 0.015 -0.064 0.072 -0.027 -0.120 -0.36 (1) -1.42 (1) -1.398 

IS 0.90 -0.024 0.053 -0.020 -0.006 0.01 0.91 0.929 

IV -3.93 -0.473 0.078 -0.127 0.038 -0.014 -0.227 -0.73 -4.66 -4.654 

8 Li(2+;l) [431] 1.16 (2) 0.226 -0.038 0.045 0.056 -0.021 0.062 0.33 (1) 1.49 (2) 1.654 

8 B(2+;1)* [431] 1.42 (2) -0.226 -0.020 -0.045 0.090 -0.033 -0.077 -0.31 (1) 1.11 (2) 1.036 

IS 1.29 -0.029 0.073 -0.027 -0.007 0.01 1.30 1.345 

IV 0.13 -0.226 0.009 -0.045 0.017 -0.006 -0.070 -0.32 -0.19 -0.309 

9 Li(§~;§) [432] 2.66 (3) 0.403 -0.076 0.141 0.108 -0.016 0.141 0.70 (2) 3.36 (4) 3.437 

9 C(|";|) [432] -0.75 (3) -0.372 0.039 -0.135 0.058 -0.031 -0.156 -0.60 (3) -1.35 (4) -1.391 

IS 0.96 0.015 -0.019 0.003 0.083 -0.024 -0.008 0.05 1.01 1.023 

IV -1.14 -0.258 0.038 -0.092 -0.017 -0.005 -0.099 -0.43 -1.57 -1.609 

9 Be(§ + ;±) [441] -1.18 (1) -0.084 0.019 -0.037 0.041 -0.009 -0.038 -0.11 (1) -1.29 (1) -1.178 

9 B(| + ;|)* [441] 2.97 (1) 0.084 -0.057 0.037 -0.005 -0.004 0.032 0.09 (1) 3.06 (1) n.a. 

IS 0.89 -0.019 0.018 -0.007 -0.003 -0.01 0.88 

IV 4.15 0.168 -0.076 0.074 -0.046 0.005 0.070 0.20 4.35 



than those obtained in the mirror symmetry based pic- 
ture. However, the error associated with this last cal- 
culation does not allow us conclude with certainty that 
implementing a CSB picture would resolve the aforemen- 
tioned 'anomaly' associated with the A = 9, T = 3/2 nu- 
clei. Within the statistics available at present, we can ar- 
gue that, in order to reproduce the experimental < a z > 
value, one would have to combine both the effect of CSB 
and the correction due to the MEC currents. As briefly 
mentioned before, we investigated the role of mirror sym- 
metry breaking also in the A = 7, 8 nuclei and found that 
the calculated m.m.'s are not sensitive to the different nu- 
clear wave functions. A trivial argument to justify this 
result is to be found in the isospin of the investigated 
nuclei: intuitively, one would expect the effect of mir- 
ror symmetry breaking to increase as AT Z = 2 T of the 



mirror nuclei becomes larger. 



B. Electromagnetic Transitions in A=6— 9 Nuclei 

In Table ITXl we report the results for the Ml transition 
matrix elements and the transition widths in nuclei with 
mass number A < 9. For these calculations — obtained 
with the xEFT EM current operator — we report only 
the IA and the MEC contributions. Exp erimental data 
in the table are taken from Refs. [3ol l36l|. The widths V 
in units of MeV are calculated as 

/ AE \ 3 

r M1 = 0.890 — B(Ml) , (26) 
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TABLE VII: Magnetic moments in nuclear magnetons of the 
A = 9, T = 3/2 systems. The nuclear wave functions are de- 
rived within a CSB picture (first block), and in a CSC frame- 
work (last two blocks). See text for further explanations. 



Nucleus(J' r ;r) 


IA 


MEC 


Total 


Expt. 


9 Li(l";|) 

3 ~ . 3 \ 
*-A2 >2> 

IS 
IV 


2.66(3) 
-0.75(3) 

0.96 
-1.14 


0.70(2) 
-0.60(3) 
0.05 

-0.43 


3.36(4) 
-1.35(4) 

1.01 
-1.57 


3.437 
-1.391 

1.023 
-1.609 


9 Li(r;§) 

9 C(§-;§)* 

IS 

IV 


2.66(3) 
-0.82(4) 

0.93 
-1.16 


0.70(2) 
-0.68(2) 
0.01 

-0.46 


3.36(4) 
-1.50(4) 

0.94 
-1.62 


3.437 
-1.391 

1.023 
-1.609 


9 Li(r;D* 
9 c(F;§) 

IS 
IV 


2.64(4) 
-0.75(3) 

0.94 
-1.13 


0.62(3) 
-0.60(3) 

0.01 
-0.40 


3.26(5) 
-1.35(4) 

0.95 
-1.53 


3.437 
-1.391 

1.023 
-1.609 



TABLE VIII: Spin expectation value (a z ) for the mirror nuclei 
9 Li and 9 C evaluated within CSB and CSC frameworks. (See 
text for explanation). 



Symmetry IA Total Expt. 



CSB 


9 Li(§" 


.3\ 9fV3- 
>2>> 


2 ) 

'2/ 


1.11 (11) 1.37 (15) 1.44 


CSC 


9 Li(|- 


.3\ 9 n /3~ 
'2>> "-^2 


.3\* 
>2> 


0.95 (1) 1.00 (1) 


CSC 


9 Li(|- 


> 2 J ' ^2 


>2> 


1.00 (1) 1.05 (1) 



where Tic is in units of MeV fm, and B(M 1) is the squared 
reduced matrix element of the m.m. operator between 
the initial and the final nuclear state divided by (2Ji + 
1), with Ji the initial state angular momentum. In the 
equation above, AE is the energy difference between the 
final and the initial state (in units of MeV) for which we 
take the experimental values as given in Refs. (35l.[36j. 

These calculations are obtained, as before, by propa- 
gating up to r = 0.8 MeV -1 with an evaluation after 
every 40 propagation steps, i.e., at intervals of r = 0.02 
MeV -1 . The analysis of the IA and the MEC contribu- 
tions are performed separately in the same fashion that 
has been implemented for the m.m.'s. 

The predictions for the A = 6, 7 nuclei as well as those 
for the A = 8, (1 + — s> 2 + ) transitions are in very good 
agreement with the experimental data. In all these cases 
the MEC corrections are needed to bring the theory in 
agreement with the experimental data. Results for the 
(3 + — > 2 + ) transitions in the A = 8 systems underpre- 
dict the experimental data, however the latter have large 
experimental errors, and thus it is difficult to reach any 



robust conclusions as to the actual level of agreement 
between theory and experiment. The transition in 9 Be 
is known with good accuracy, but the predicted width 
is lower than the experimental data although the error 
bars almost touch. We also report a prediction for the 
(h — > | ) Ml transition in 9 Li which has not been mea- 
sured yet. We did not calculate the 9 C transition to its 
unbound i state. 

The magnetic transition densities in IA as obtained 
from the VMC starting wave functions are shown in 
Fig. |U As before, the red upward-pointing triangles 
are the contribution from the proton spin term, the blue 
downward-pointing triangles are from the neutron spin, 
the green diamonds are from the proton orbital term, and 
the black circles are the total IA contribution. For the 
lithium isotopes, the transitions are predominantly from 
the proton spin term, i.e., these are almost pure proton 
spin-flip transitions. For 7 Be and 8 B, the neutron spin 
term is the most important, but with some contribution 
from the proton spin and orbital terms. The neutron 
spin-flip is also the biggest term in the 9 Be transition, 
but here the proton orbital piece is almost the same size. 

Finally in Table [Xj we show IA results for the electric 
quadrupole matrix elements and the associated transition 
widths. The latter in units of MeV is 

(AE \ 5 
— J B(E2) , (27) 

where B{E2) is the square of the reduced matrix element 
of the electric quadrupole operator given by 

p 1A = J2^r^Y 2 (r l ) , (28) 

i 

where Yj, is the spherical harmonic. The I A picture pro- 
vides a decent description of the two experimental data 
points that are available, which might possibly be im- 
proved by the inclusion of two-body effects. This topic 
has not been addressed in this work although effort in 
this direction is underway. 

The results discussed in this section are summarized 
in Fig. [7] for EM transitions whose widths are known ex- 
perimentally. We observe that in most cases, the agree- 
ment with the experimental data is excellent, and that 
the MEC contributions are crucial for the B(M1) cases. 



V. CONCLUSIONS 

In this work we have reported GFMC results for EM 
moments and transitions of A < 9 nuclei. The calcula- 
tions of m.m.'s and Ml transitions account for the effect 
of two-body EM currents, for which we considered two 
models: i) the SNPA model described in Refs. |2 l . [lCI. and 
ii) the pionful ^EFT EM operator derived in [l]]-[l3( . The 
goals of this work were to continue the study initiated in 
Refs. dHH], by extending the calculations to systems with 
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TABLE IX: Matrix elements in units of nuclear magnetons and widths of Ml transitions in A = 6-9 nuclei which account for 
the xEFT current operator up to N3LO. IA and MEC contributions are shown. 



(J? H 


■ JJ) 


Ml and T IA MEC Total 


Expt. 


6 Li(0+ 


-+1+) 


Ml 3.63 (1) 0.38 4.01 (1) 
T (eV) 6.90 (2) 8.41 (3) 8.19 (17) 


7 Li(I- 


->■ D 


Ml 2.66 (1) 0.47 (1) 3.13 (2) 
T (10~ 3 eV) 4.47 (5) 6.19 (8) 6.30 (31) 


7 Be(r 


. 3-' 
2 - 


1 Ml 2.31 (2) 0.41 (1) 2.72 (2) 
T (10~ 3 eV) 2.44 (4) 3.39 (6) 3.43 (45) 


8 Li(l+ 


^2+) 


Ml 3.47 (4) 0.74 (2) 4.21 (5) 
T (KT 2 eV) 4.4 (1) 6.5 (2) 


5.5 (1.8) 


8 B(1+ - 


^2+) 


Ml 3.17 (5) 0.67 (2) 3.84 (6) 
T (KT 2 eV) 1.8 (1) 2.6 (1) 


2.52 (11) 


8 Li(3+ 


^2+) 


Ml 0.98 (6) 0.20 (5) 1.17 (8) 
T (10~ 2 eV) 1.8 (2) 2.6 (3) 


7.0 (3.0) 


8 B(3+ - 


-> 2+) 


Ml 1.31 (6) 0.23 (5) 1.56 (8) 
T (KT 2 eV) 3.5 (3) 4.9 (5) 


10 (5) 


9 Li(r 


-+ D 


Ml 2.28 (3) 0.36 (4) 2.64 (5) 
T (10" 1 eV) 5.9 (2) 7.9 (3) 


n.a. 


9 Be(|- 




1 Ml 1.42 (3) 0.20 (2) 1.62 (4) 
T (10~ 2 eV) 5.6 (3) 7.2 (4) 


8.9 (1.0) 



A > 7, and to test the x^FT two-body EM current oper- 
ator within a hybrid context. Both models describe the 
long-range behavior of the two-body EM current in terms 
of OPE contributions. These pseudoscalar terms consti- 
tute the major contribution to the total MEC correction. 
The models also include, albeit in different formulations, 
the effects due to currents involving A isobar degrees of 
freedom. While the SNPA current does not involve free 
parameters, the %EFT EM operator invokes a number of 
unknown LECs which have been fixed to reproduce EM 
observablcs in the A = 2,3 nuclei. This additional free- 
dom is probably responsible for the closer agreement with 
experiment given by the %EFT formulation compared to 
the SNPA model. In particular, both the isoscalar and 
isovector MEC corrections for the m.m.'s are closer to 
experimental data when calculated with the xEFT EM 
currents. Nevertheless, we find that the two models are 
in good qualitative agreement and both support the ne- 
cessity of adding MEC corrections to reach agreement 
with the experimental data. 

In view of the study presented in Ref. [45| and ref- 



erences therein, we have paid special attention to the 
9 Li, 9 C mirror nuclei. The experimental value for the 
isoscalar spin expectation value, (a z ) = 1.44, has been 
considered 'anomalous' and various explanations, includ- 
ing a broken mirror symmetry, have been suggested. We 
find that calculations of (a z ) obtained assuming mirror 
symmetry are close to the single-particle estimate of 1, 
even if MEC contributions are included. When the mir- 
ror nuclei wave functions are constructed individually, in- 
cluding the appropriate Coulomb differences in both the 
starting VMC trial function and GFMC propagator, and 
thus breaking mirror symmetry, both the IA and MEC 
components of (a z ) are increased in the right direction, 
giving a result that is consistent with the experimental 
value, and indicating that both broken symmetry and 
MEC contributions arc required. However, the error as- 
sociated with this last calculation of the (a z ) value is too 
large to make a definitive statement at this stage. 

Finally, we have studied a number of EM transitions 
induced by the Ml and E2 operators. After including the 
MEC contributions, the calculated Ml transition widths 
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FIG. 6: (Color online) Ml transition density in nuclear magnetons per mi 3 for selected nuclei (see text for explanation). 



are in excellent agreement with the experimental data for 
the A = 6, 7 nuclei and the A = 8 (1 + — >2 + ) cases, while 
theory undcrpredicts the data in the A = 8 (3 + — > 2 + ) 
and A = 9 cases. However, the latter A = 8 transi- 
tions have large experimental errors, so new precise pre- 
cise measurements for these as well as for the as yet un- 
measured 9 Li(| — > | ) transition would be most useful. 
For the E2 transitions we provide only I A results. While 
two-body corrections to the IA charge operator have been 
derived in both SNPA [Sill and %EFT H3, neverthe- 
less their leading contribution due to OPE is expected 
to be small, because the associated operator vanishes in 
the static limit. The present IA calculations appear to 
describe satisfactorily the E2 transition widths for both 
the 7 Li(±~^§~) and 9 Be(§~-»§~) cases. 
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